A General Approach to Casimir Force Problems Based on Local Reflection 

Amplitudes and Huygen's Principle. 
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In this paper we describe an approach to Casimir Force problems that is ultimately generalizable 
to all fields, boundary conditions, and cavity geometries. This approach utilizes locally defined 
reflection amplitudes to express the energy per unit area of any Casimir interaction. To demonstrate 
this approach we solve a number of Casimir Force problems including the case of uniaxial boundary 
conditions in a parallel-plate cavity. 

I. INTRODUCTION 

Since Casimir's pioneering work in 1948 many advancements have taken place in the field of Casimir force theory 
and its experimental confirmation. Casimir's original work described the force that operates between uncharged 
parallel conductors arising from the vacuum energy of the electromagnetic field. The theory has since been extended 
to isotropic dielectrics [l|-|3| and anisotropic materials [J, Q. Experimentation in the field began in the 1950's with the 
efforts of Overbeek and Sparnaay[g and of Derjaguin and AbrikosovafT] and has continued to the present day with 
ever increasing precision [8|-[ll|, covering a wide range of separations from lOnm— 1/xm. 

However, 63 years after the birth of this field of study, many problems still remain. These problems include the 
modelling of material properties, the effects of temperature, and the approach to complex geometries among others. 
In addition to these out-standing problems, there exists another problem: the lack of a common approach. From 
mode counting methods to Green's functions to stress-energy tensor formalisms, each constitutes a unique method 
for calculating the Casimir force using exclusive nomenclature leading to unique interpretations of the phenomenon. 
Although most of these approaches yield the same or similiar predictions, it is often difficult to know exactly how to 
use any particular formalism to solve a new problem or assess the domain of its applicability. In an effort to redress 
this issue a number of excellent books have been written that detail many of these approaches [12 - [l4| . Herein we 
develop an approach that covers a wide-range of Casimir problems and we hope to describe it in detail so as to be 
easily adapted to new problems. 

In this paper we will calculate the Casimir force using a scattering approach and subsequent mode counting method. 
Our basic formalism is not specific to any geometry or boundary conditions; we simply assume that we have two 
interacting boundaries as shown below and note that the right-bound waves in the cavity are generated by reflections 
of all the left-bound waves off the surface 5*1 and vice-versa. The amplitude of the right-bound wave and left-bound 
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FIG. 1. Two interacting surfaces Si and S2 of arbitrary geometry. The field between the surfaces will be denoted as the sum of a 
right-bound wave and a left-bound wave. 



wave at points xi and X2 on surfaces Si and 5*2 respectively are related by a generalized Huygen's Principle, namely, 
that the amplitude of any polarization state at a particular location is given by the sum of all incident wavelets added 
up in a phase sensitive manner: 



^fl(x2) = / .g(x2,xi)r(x2,xi)i?i*L(xi)dS'i 

J Si 
*l(xi) = / g(xi,X2)r(xi,X2)i?2*i?(x2)dS'2 

Js^. 



(1) 



(2) 



Here i?i and R2 are n x n reflection matrices with dimension equal to the number of polarization states of the field 
under consideration (for the electro-magnetic field these are 2x2 matrices) and gT equals the gradient of the Green's 
function. Specifically, 

where F, hereafter referred to as the transport functions or transport matrix, encodes only the phase changes and g 
encodes only the modulus changes suffered by the wavelet during transport from some point x to some other point 
x'. Both g and T are n x n matrices. The properties of Green's functions and Eq. |3] requires g to satisfy: 



g(x,x')5(x',x")d5d5"=/ (4) 

Substituting equation 2 into equation 1 we get, 

*fl,(x2)=/ / g(x2,xi)r(x2,xi)i?ig(xi,x'2)r(xi,x'2)i?2*fl,(4)d^id^2 

J Si J S'^ 

Noting that gVRgTR is invariant under cyclic permutation we can rewrite the above as: 

*i?,(x2) = / / g(x2,Xi)5(xi,x'2)r(x2,xi)i?ir(xi,x'2)i?2*i?.(x2)(iS'id5'2 

J Si J S' 



(5) 



(6) 



r(x2,xi)i?ir(xi,x2)i?2 = | 5 (7) 

The operator on the left-hand side of Equation [7] can take on either of two values, or I. If the operator equals 0, 
^ r{'^2) = for all X2. Therefore, in order for a non-trivial solution to exist, 

det(I-r(x2,xi)Rir(xi,X2)R2)=0 (8) 

An eigenmode of the cavity exists for any choice of oj, k, xi, or X2 which can satisfy equation [8l 

In order to calculate the Casimir interaction energy one must sum over the energies of all such eigenmodes, a task 
most easily performed by utilizing the argument principle: 
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(9) 



In the above expression A = det(I — r(x2,Xi)Rir(xi,X2)R2). The transport functions, F, have elements proportional 
to e"'*''', making the above integral oscillatory. To deal with this we can perform a change of variables, ^c = —iuj/c. 
All the elements of F will then be proportional to exponentially decaying functions and consequently the integral over 
the semicircle C4. goes to zero. After integrating by parts the expression for the energy carried by all the modes in 
the cavity is given by: 

+■ />C30 

E=-^Y. ln(det(I-F(x2,Xi)RiF(xi,X2)R2))dCc (10) 

paths 

where the sum is over all unobstructed straight line paths connecting the two surfaces. 

II. COMPARISON WITH PREVIOUS RESULTS 

In this section we will demonstrate how to use this approach by solving a few specific cases and showing that this 
approach yields results consistent with past works. We will start with the simplest case, the Casimir case described 
by a parallel-plate cavity of infinite lateral extent L^ ^ Ly ^ 00 and finite spacing L^ = a, the walls of which are 
perfect conductors, Ri = R2 = I. The cavity itself is empty with e = 1, a fact that allows us to say the transport 
functions are diagonal (no polarization mixing upon transport). Moreover, the transport of an eigenmode across the 
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cavity is invariant under the change of polarization, this aUows us to write the transport functions as: Fr = e*''''!, 
where r = (r^., Ty, a) is the vector between xi and X2. From this we can write equation 1101 as: 

E^~— f dk f dr f ln(det(I - e-"'"-e-*'(-'-)l))d^e (H) 



Here k' = (fcj_ cos0, k± sincj), —i\JQ + A:^) and k = {k± coscj), k± sm4>, i \JQ. + fc^. ) ^^^ ^^^ ^^^^ ^^^ right bound wave 
vectors, respectively and the dk and dr integrals are a symbolic representation of the sum over all the connecting 
paths. We can express the integrand in a different form by noting that for any 2x2 matrix. A, the following is true: 

det(I- A) = l-tr(A)+det(A) (12) 

Using equation 1121 and the definitions of the wave- vectors we can rewrite equation [TT] as: 

he f°° k\ dk\ 



E = -— ^i^^ / dr / ln(l - 2e-2-V?H^ + e-^'^V^?+^)d^, (13) 

47r Jo 27r ' 



— 00 



E/A = -^j k^dk^j ln(l-e-2«V«I^)dC. 



(14) 



Here we have written E / A as the energy per unit area in order to eliminate the integral over dr which symbolically 
represents an area integration the details of which will be addressed in the next section of this paper. We can directly 
evaluate the above equation by expressing the natural log as a series and performing the integration term by term: 

/•oo /'Oo 00 



E/A=^ k^dkj d^^X^V^anVlTIT (15) 

27r2 Jo Jo t^n 



Tjg °° r'^/2 I'ca -2anp 

= ---^limV/ / p^sinOdedp (16) 

2n2s^o^Jo Js n ' "^ 

here we have performed a change of variables: ^c — pcos0 , k± — psuiO. The integration over 9 can be carried out 
quickly to give a value of 1. The lower limit of integration of p has been written as S to eliminate any controversy over 
expanding the natural log around a discontinuity. Integrating over p by parts and taking the limit as (5 — > yields 
l/{4a^n^), so that: 



E/A = ^V^ (17) 



87r2a3 

n— 1 

The sum in the above expression is well known and has a definite value of C(4) = 7r'*/90. This yields the value for the 
energy per unit area: 

i^...(a)^-^ (18) 

in conformity with the standard Casimir result [l2|. 

III. PARALLEL-PLATE CAVITY OF FINITE SIZE 

In the last section we calculated the Casimir energy per unit area of a parallel-plate cavity of infinite extent. 
Although this is an interesting result it is only applicable under the condition that the separation of the plates is 
much less than the size of the plates. Below we calculate the Casimir energy per unit area of a cavity of finite size. 

For simplicity let us once again consider the case of a perfectly conducting parallel-plate cavity. This time however 
the side-walls will only occupy the region (—6, b) in the x- and y-directions. The separation of the plates will once 
again be a as shown in Figure [21 To solve this problem lets start with equation [TT] and write the integral over r as an 
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FIG. 2. Two parallel plates separated by a distance a in the z-direction. The plates are squares of side-length 26. The plates are aligned 
such that the bottom left corners of both plates are given by the xy-coordinates (— b, —b) and the top right corners are given by (6, b). 

integral over either surface weighted by the probabihty that the set of isotropicahy distributed path vectors emanating 
from each differential area element is intercepted by the other surface. 
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/ k^dkj_ / ln(l - e-2V«c+fcl)^^^ 
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(19) 



The factor of tt in the denominator of the r integration is the normalization factor that is derived by noting that for 
a plate of unit area placed parallel to an infinite plane the spatial part of the integral should be unity, in order that 
Eq. [12] agree with Eq. [TH]in the limit & — > oo. Writing Eq. \W\ explicitly: 



E = - 



TT^hc 1 
720a3 ^ 



b/a rb/a rb/a rb/a 
-b/a J -b/a J -b/a J -b/a 



dxidyidx2dy2 



{x2 - xif + (y2 - yif + 1 



(20) 



here we have performed the integration over the wave properties in the same manner as the preceding section and 
scaled the integration variables by a. From this we can see that Eq. [20] differs from Eq. [18] solely in the manner in 
which the area is calculated. To demonstrate the effects of finite size on the overall Casimir energy we have plotted 
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Plate Separation {ntn) 



FIG. 3. Log-Log graph of energy vs separation. This graph was generated for square plates of side-length 10 fim separated by distances 
between .1 fira and 100 ^m. The dashed line represents Casiinir's expression Eq. llSl mulitplied by the area of Plate 1. The solid line 
represents our expression Eq. 1201 



the Casimir energy vs separation for square plates of side 26 = 10/im over separations of .1 ^m to 100 /im. The 
dashed line in the graph was calculated using Casimir's expression Eq. 1181 simply multiplied by the area of Plate 1, 
the solid line was calculated using our expression Eq. ^U\ From Figure [3] it is apparent that for cavities with lateral 
dimension much larger than their separation (b/a >> 1) our expression approaches Casimir's expression, however for 
separations larger than the side- length (b/a < 1) our expression begins to drop off quicker than Casimir's expression. 
This behaviour is to be expected, as the plates should begin to appear as points rather than planes for large separation, 
leading to an additional 1/a^ behaviour. 



IV. PARALLEL-PLATE CAVITY DESCRIBED BY UNIAXIAL BOUNDARY CONDITIONS 

An interesting case is that of a parallel-plate cavity with uniaxial boundary conditions as illustrated in Figure |4l 
Such a cavity would have side- walls made of a uniaxial material, for example a wire grid polarizer or a graphite crystal 
cut with in-plane optical anisotropy. When two such plates are placed close together with their optic axes rotated by 
an angle x with respect to each other, both an orientationally dependent normal force and a torque tending to align 
the optic axes of the plates have been predicted Q [1^1 [2J| . In order to address this Casimir problem in the formalism 

i?.[Y] 
e2± / 
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FIG. 4. Two parallel plates separated by a distance a in the z-direction. The coordinate axes are oriented to coincide with the principle 
axes of Plate 1 (bottom plate) with the x-direction coinciding with its optic axis. Plate 2 (top plate) is rotated about the z-axis by an 
angle x relative to the Plate 1. The dielectric tensors, e, associated with each plate are also given. 

presented in this paper, let us start with equation 1101 relevant to a parallel-plate cavity of infinite extent: 

E/A ^-^ dk ln(l - tr(RiR2)e-2V«?+kl + det(RiR2)e-^V«?+'^i)dCc (21) 

2vr J Jq 

We need only to define Ri and i?2 and proceed as before. In order to multiply them is this fashion these refiection 
matrices need to be in the same basis. We will write our reflection matrices in the TE-TM basis. Fortunately for all 
parallel-plate cavities the basis vectors that define the TE-TM basis for a given wave-vector on one surface are the 
same as for the other surface. Thus i?i and i?2 will automatically be in the same basis. Were this not the case we 
would write Ri = Bq^^lTZiBl^^q, where TZ refers to the reflection matrix in a local coordinate basis. L and G refer 
to a local and global basis and B indicates a change of basis matrix which could be a function of position. 

The reflection matrices are 2x2 matrices the elements of which describe transition amplitudes between the polar- 
ization states of the incident and reflected waves. 

\ riEM riMM J \ r2EM r2MM J 

Where Ri and i?2 represent the reflection matrices of Plates 1 and 2 respectively. The functional definition of these 
matrices is to describe the reflected wave in terms of the incident wave. So, if we let -0^ = (1 0) represent an incident 
transverse electric wave with unit amplitude, we can write the reflected wave as V'r = Ri'4^i = [tiee fiEM)- 

In order to obtain an explicit expression for the matrices we must find the general solution for a wave in each of 
the three regions (the interior of the cavity and inside plates 1 and 2). We then have to match up these solutions 
with Maxwell's equations providing the boundary conditions. In the region between the plates (Region 3) the general 
solution of Maxwell's equations is a linear combination of transverse plane waves described by an angular frequency 
w, wave vector k, and polarization unit vector A. For any choice of uj we can write the solution as the sum of forward 
and backward travelling waves, each of which have two independent polarizations. Thus the solution for any of w will 
be the sum of four waves, written in terms of the vector potential as: 

A3 = (aaAs -I- a4A+Af )e'('^-"-"*' + (asA^ + ae\_M)e'^'''-'"^'^ (23) 

In the above equation a^ is the amplitude of the z*^ wave. The polarization and wave vectors can be written explicitly 
as: 

k± cos{(j)) \ / k± cos(0) \ / sin(0) \ / cos{(l))p3 

k±_ sin((/)) , k' = /ci sin((/)) , A^ = - coii{<j)) , A±m = , ^ ^, sin(0)p3 | (24) 



In the above equations /c^ is the magnitude of the component of k perpendicular to the z-axis and p^ ~ ^/eaCc + ^1 
where we have performed a scaled Wick rotation as before, ^c = —iijj/c. 

The general solutions in the interior of Plates 1 and 2 (Regions 1 and 2, respectively) [16] are similar but given in 
terms of a linear combination of ordinary and extraordinary waves denoted by the subscripts o and e. In contrast 
with Region 3, the solutions in Region 1 and 2 will have two waves each instead of four because we assume Regions 
1 and 2 to extend to infinity thereby eliminating any inbound waves. The solutions in Region 1 and 2 are then: 

1 — aiAioC ^ '-" ' + a2Aiee ^ ' , A2 — ajA2oe ^ -° ' + asA2ee ^ " ' 1,/oj 

It's well known that the tangential components of the wave vector are left unchanged by refraction and as such, only 
the normal component of the wave vectors and the three components of the polarization vectors in the above expression 
need to be worked out explicitly. This can be done by solving the set of four coupled algebraic equations; (1): the 
norm of the polarization vector equals one (A^ — 1), (2): the transversality of the wave (k-D = 0), (3): the definitions 



of the ordinary (D • 6 = 0) and extra-ordinary waves(D • (kg x 6) = 0), (4): Snell's law I k^ — ±i\/S,c (e^) 
Where D = eE is the electric displacement. The resulting wave and polarization vectors can be written explicitljt^: 




kioz = -ipi = -iJe^^Q + fci fcie. - -i~Pi = -^\hnG + ^i + ^1 7^ - 1 cos2[0] (26) 



= iP2 = i^Je2±ec + fci fc2e. -*P2 = ye2||e? + A:i + fcl(f^-l)cos2[0 + x] (27) 

\ ^ /eMklsin^[cl^]-pl)\ ^ ( sin[x]p2 \ 

Aio = ( -Pi Ale = fc5^ei||cos[0]sin[0] A20 == cos[x]p2 (28) 

ik^s\n[4)\ j y ifc^eijl cos[0]/5i j y ifc_L sin[(/) + x] / 

/ k\ sin[0 + x](e2|| cos[(/) + x] sin[x] - e2± cos[x] sin[0 + x]) + e2_L cos[x]pi \ 
Ase = - k\ sin[0 + x](e2|| cos[0 + x] cos[x] + £2^ sin[0 + x] sin[x]) - e2_L sin[x]pl (29) 

\ ik^e2\\cos[(j) + x]p2 ) 

We can now use the fact that that the tangential components of the electric and magnetic fields are continuous across 
the boundary to solve for the reflection coefficients. The matrix elements for Ri are given explicitly in Appendix A. 
From the reflection coefficients it is obvious that the Casimir force depends on the angle of orientation, x, which can 
not be decoupled from the azmuthal angle (j). Thus the integration over dk in equation [2l1 should be replaced with 
integrals explicitly over dk±^ and dcj): 

E/A = ^/ k^dk^ d(t> ln(l - tr(RiR2)e-2aV«I+kl + dct(RiR2)e-4"V?I+kl)dee (30) 

Two cases of the above expression are of particular interest, the case where e^n = ei±_ — ti (isotropic case) and the 
case where £3 = ein = 6211 = 1 and ei^ = e2i_ = 00 (perfectly anisotropic case). We consider these cases below. 

A. Isotropic Case 

It can be shown that in the isotropic limit the reflection matrices given in Appendix A reduce to the isotropic 
reflection coefficients. In this case the reflection matrices are: 

C -P1+P3 \ / -P2+P3 Q \ 

P1+P3 i?9 = ''2+P3 

n eaPi-eiPa I i ^"-2 I q £3P2-e2P3 I 

e3Pl+ciP3 / \ e3P2+C2P3 / 

the elements of which equal the standard reflection coefficients for the isotropic case. The integrand in equation 1211 
thus reduces to: 

ln[l - {riEr2E + riAfr2Af )e^2p3a ^ {riEriMr2Er2M)e'^''"'] (31) 



^ Note that the normahzations on the polarization vectors have been suppressed. 



Here we have dropped the second letter in the subscripts as the matrices are now diagonal. If we let ei = £2, 
{riE = r2E, Tim — '''2m), Eq. 1301 simplifies further to yield the standard Lifshitz expression: 

E/A = -—j j k^H{l-rle-^^P'^){l-rl,e-^^P-)]d^,dk^ (32) 

B. Perfectly Anisotropic Case 

Another interesting limit of our expression is the totally anisotropic case. In this limit we let e^ — ^i\\ — ^2|| ==1 
and we will let eix and e2X go to infinity. The resulting reflection coefficients for Plate I are then: 





k]_sin^ <p + ^l^ 


(1 -2cos^0) - 


-^fei+C2^feisin2 + 42 




k^ sin^ 


+ c? + y^I 


+ ??^fclsinV + ?? 




aCc^fci 


+ ?2 cos[0] sin[0] 




fci 


1(1 -2cos^0) 


/^l + CcV^i "'"'•* + «' 




-^/kl+^^^^klsin^4> + e. 



kl sin2 + e? + ^fe^ + C? ^/ci sin^ + e? 



(33) 



(34) 



(35) 



The reflection matrix for Plate 2 is identical to that of Plate I with the replacement oi (/)—>■ (j) + x- As the interaction 
energy depends on the product of these matrices, it is demonstrative to see the form of the product under a few simple 
cases, namely under parallel (x — 0) and perpendicular (x — 7r/2) alignment. 

10 

^1-^21x^0,0^0= I „ f ^c-\/k^_L+ic\ ) ^1^21x^0,0^71/2= ( n 1 ) ^"^^^ 



\ / 

^1^2|x-^7r/2,0^O = ( Q -Cc + x/fcj+Cg , -Rl-R2|x^7r/2,0^7r/2 = I q -gc + x/fcj+gg | (37) 

i. + y/kl+il J \ i. + y/kl+Ca 

From this it is easy to see that for perpendicular alignment the trace of the product of the reflection matrices is 
reduced with respect to parallel alignment, however it does not go zero for all cases. Thus, one would intuitively 
predict that the interaction energy will not vanish for perpendicular alignment. In order to obtain the energy per 
unit area as a function of the orientation angle (x) we have evaluated eq. [301 using the reflection matrices from eqs. 
[551155] and plotted the values in Figure [5] for a separation of 1 micron. The values have been represented in units of 
the energy in the Casimir case (Eq) and fit to a cos^[x] function. 

The 180° symmetry in the orientational dependence could be of particular interest to experimentalists in the field. 
An experiment could be constructed whereby one of the plates would be rotated relative to the other at a set frequency. 
The force between the two plates would thus vary at twice the frequency, while all other signals should remain constant 
or vary at the rotation frequency. Using this method the Casimir force should be easier to isolate from other signals 
infiltrating such an experiment. The dependence of the energy on the angle of orientation will also produce a torque: 
M = ^. This torque will exhibit the same 2x dependence as the normal force, however it will be maximum at 
X = 7r/4 instead of x = 0- Such a torque would produce a displacement orthogonal to that of the normal force and 
could in principle be used as another method of measuring the macroscopic effects of the quantum vacuum. 

V. FINAL REMARKS 

In this paper we have developed a general approach to Casimir problems and have used this approach to solve a 
number of examples. We have shown that this approach yields Casimir's and Lifshitz' expressions in the appropriate 
limits. Using this approach we have also calculated the Casimir force between anisotropic conductors and shown that 
this leads to the predicted orientational dependence. It is our hope that we have described this approach in enough 
detail so that others will find it easy to adopt the formalism in their own calculations. 

The authors wish to thank Girish Agarwall for his help and insight. This work has been partially funded by an 
NSF Grant #61018 allocated to R.C. 
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FIG. 5. The Normalized Casimir Energy per Unit Area as a function of orientation for a cavity formed with two perfectly anisotropic 
plates. The values of the energy have been scaled with respect to the (isotropic) perfectly conducting case. The bluc(top) line represents 
the fit, fix) = 0.2286 + (0.5235 - 0.2286) cos^fx]. 



Appendix A: Reflection Matrices 



This formalism requires explicitly defined local reflection coefficients, in order to obtain them for boundaries defined 
by in-plane optical anisotropy we will match the general solutions given in Eqs. 1231 and [^ using Maxwell's equations 
as the boundary conditions, which state that the tangential components of the electric and magnetic fields must be 
continuous across the boundary. Explicitly, the boundary conditions for Plate 1 are: 



ijI,E^Ex + ^I,-M^-Mx + i^RE^Ex + '0_R, + A/ A+Ma; — ^To^ox — i^Te^ex = 



(Al) 



fpI,E^Ey + i'l.-M^-My + i^RE^Ey + i^R^ + M^+My — 4'To^oy — i'Te^ey = 



(A2) 



'(/'/,£(kx A£;):j + ?/'/,-A-/(kx A_M)a;+'0i?,B(k' X Xe)x +1pR, + Mi'ii-' X X+m)x -IpToiMlo X Xo)x -IpTeO^le X Ae)^: = (A3) 



V'/,B(k X XE)y + 1pI,-M{^ X A_Af )y +^fiB(k' X Ab)^^ + ^fl,+M (k' X X+M)y - IpToCi^-lo X Xo)y - ^Te(kle X Xe)y = (A4) 

here ipij is the amplitude of the ith wave (i standing for Incident, Reflected or Transmitted) with jth polarization 
(j standing for transverse Electric, transverse Magnetic, ordinary, or extraordinary) and Xjk is the fcth component 
(fc standing for the x or y directions) of the jth polarization unit vector. From these four equations, the amplitudes 
of the reflected waves can be written in terms of the incident waves, convienently arranged in a 2 x 2 matrix whose 
elements are given explicitly below. 

riEE = [fci sin2[(/i]p3 (p3 - pi) {A + ei±pl} + plpsPi {A + ei±p?} + pi {fc^^ei^ sin2[</)]/93 {pi - k^) 
~A{kl sin2[0] + cos[20]p§)pi + ei±p3 [kl - pfj p\ - ei^ (k\ sin2[0] + cos[20]p§) pf] 



'A{k\ sin^[0] + cos[20]p§)pi + eu^pa {k\ - p§) pi - ei±_ {k\ sin^ 
-pi {k\ cos^[(J3\ (fci(ei|| - eiA_)sva^[(J3\ + ei^Pi) 



-pf {/eicos"[0J (fel(ei|| -ei±)sm"[0J + eixplj 

+pI {k\ siT?[4>] (eii cos[24>] ~ 2e^ cos2[0]) + ei^ cos[2(/)]p2) }] ±- 



(A5) 



riEM = ^ cos[0] sin[0]p3^p2 - k^ [k\{ei\\ - €i^)sm^[(j)\ {k^ -pi) + Apipi + ei^ {k^ - pf) pi + ei±pipf\ (A6) 



riME = 1^ cos[0] sin[0]pip3^p2-A:2^(pi - pi) [A + ei±p?] 



(A7) 



riMM = [k\s\T?[4)\{-Apl + ei^pip'i (p§ - k\) + pj (fci(ei|| - eix) cos2[0] - (2ei|| cos2[0] - ei^ cos[2(/)]) p§)} 
+A {/ci sin2[(/>]pi - {kl sin2[(/)] - pf) ^3 + cos[2(/)]pip^} pi 
+ei± {-fci sm^[(j>]pl + pip3 {kl - pI) + pI {kl cos2[0] - cos[2(/.]p§) } pi 



+ei±_ {kl svc?[4>]pi + {pI - kl sin2[0]) pg + cos[2(/)]pip2} p3l 1 



where, 



(A8) 



A= -kl{€n -{en -eii)sin2[ 



^1 = iPi + Pa) [^i sin2[<^] (p3 + pi) {A + ei^pj} + pi [ei^pi {/c^ sin2[0] - p2} 
-P3P1 {A + ei_LP?} + A:i cos2[0] {fc2^(ei|| - ei_L) sin^[0] + ei±p?}]] 

Here riy represents the amphtude of the reflected wave with j polarization in terms of an incident wave of i polarization 
at the surface of Plate 1. The reflection coefficients for Plate 2 can obtained by the simple replacement of — >■ 4> + Xj 

Pi -^ P2, Pi -^ P2, ei|| -^ £211 , eiJ_ ^- e2_L- 
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